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SUMMARY 

Nonlinear wave-wave interactions in turbulent jets were investigated based 
on the integrated energy of each scale of motion in a cross section of the 
jet. The analysis indicates that two frequency components in the axisymmetric 
mode can interact with other background frequencies in the axisymmetric mode, 
thereby amplifying an enormous number of other frequencies. Two frequency com- 
ponents in a single helical mode cannot, by themselves, amplify other frequency 
components. But combinations of frequency components of helical and axisymmet- 
ric modes can amplify other frequencies in other helical modes. The present 
computations produce several features consistent with experimental observations 
such as (1) dependency of the interactions on the initial phase differences, 

(2) enhancement of the momentum thickness under multifrequency forcing, and 

(3) an Increase in background turbulence under forcing. In a multifrequency- 
excited jet, mixing enhancement was found to be a result of the turbulence 
enhancement rather than simply the amplification of forced wave components. 

The excitation waves pump energy from the mean flow to the turbulence, thus 
enhancing the latter. The high-frequency waves enhance the turbulence close 
to the jet exit, but the low-frequency waves are most effective further 
downstream. 


INTRODUCTION 

The study of artificially excited shear flows may provide insight into the 
physical mechanisms governing the onset and evolution of turbulence in shear 
flows. Such studies can provide a guide for formulating new predictive turbu- 
lence models and for controlling lift and drag as well as turbulent mixing and 
transport. Motivated by understanding laminar-turbulent transitions, Miksad 
(refs. 1 and 2) excited a shear layer with two frequencies of comparable growth 
rates. The two imposed waves interacted nonlinearly, generating the sum and 
difference modes and the harmonics and subharmonics of the forcing frequencies. 
Some of these components were found to grow and interact with others, generat- 
ing additional components. Finally an enormous number of waves were generated, 
leading to a complete transition to fully turbulent flow. Because of the com- 
plexity of this experiment, most later experiments concentrated on single- 
frequency excitation. 

Single-frequency excitation has now been studied by many researchers: 
for example, Binder and Favre-Marinet (refs. 3 and 4), Ho and Huang (ref. 5), 
Mankbadi and Liu (ref. 6), Mankbadi (ref. 7), Winant and Browand (ref. 8) and 
Zaman and Hussain (ref. 9), to mention only a few. An excellent review for 
the work in this area is given in an article by Ho and Huerre (ref. 10). In 
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general, these investigations indicate that the shear-layer spreading rate can 
be controlled effectively by forcing. 

The possibility of using two-frequency forcing to extend the range over 
which the flow can be controlled has been explored by Wygnanski and Petersen 
(ref. 11), Raman and Rice (ref. 12), and Bradley (ref. 13). These observations 
indicate that forcing at two frequencies makes the shear layer grow faster than 
forcing at either frequency alone at the same forcing level. Thus, combined- 
mode forcing is not only significant to the fundamental understanding of the 
laminar-turbulent transition, but is also significant to the technologically 
important problem of mixing and turbulence control. The importance of two- 
frequency forcing has motivated several studies to understand the wave-wave 
interactions in two-frequency excitations. Numerical simulations of temporal 
instability in a shear layer (refs. 14 and 15) and of spatial stability in a 
laminar axi symmetric shear layer (ref. 16) show that the wave-wave interaction 
process is highly dependent on the phase alignment. This was confirmed exper- 
imentally by Zhang, Ho, and Monkewitz (ref. 17) for the shear layer. For a 
round jet, Arbey and Ffowcs-Wi 1 1 iams (ref. 18) observed a phase-dependent sup- 
pression of the first subharmonic in the pressure spectrum when the jet was 
excited at the natural frequency and its first subharmonic. Furthermore, the 
experimental observations of Raman and Rice (ref. 12) show that for high levels 
of excitations, the phase development of the stability waves is a highly non- 
1 inear process . 

This report presents formulations for multifrequency Interactions in the 
axi symmetric or azimuthal modes. These wave-wave interactions are coupled 
with their interactions with the fine-grained background turbulence. The pur- 
pose is to provide a model for a multifrequency excited jet and to Increase 
the understanding of the later stages of the laminar-turbulence transition. 
High-amplitude excitations are emphasized so that the maximum effect on jet 
mixing can be observed. The nonlinear development of the wave amplitudes are 
coupled with the nonlinear development of the phase angles, which Raman and 
Rice's (ref. 12) observations indicate to be significant in high-amplitude 
excitations . 


ANALYSIS 

The problem considered here is that of a turbulent jet issuing from a noz- 
zle of radius R, the Mach number is small, and the compressibility effects 
are negligible. The coordinate system is (x,r,<j>) where x is the streamwise 
direction, r is the radial direction, and 4> is the azimuthal angle. The 
corresponding velocity components are u, v, and w. Velocities are normal- 
ized by the jet exit velocity U e , distances by the nozzle radius R, and the 
pressure by pl)2, where p is the fluid density. The time t is normalized 
by R/U e . e 

The_fluid motion is split into three kinds of motion: a time-averaged 

motion Ui(x); a periodic, organized, large-scale wavelike structure Uj(x,t); 
and a background, fine-scale random turbulence u!(x,r). Thus, 


u.(x,t) = U.(x) + u^(x,t) + u^(x,r) (1) 
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and the pressure, P, is similarly split. An overbar denotes the usual time 
average of the flow quantity. The periodic, wavelike component is taken to be 
represented by a finite number of frequency and azimuthal components. Thus, 


t 


5<*.t> - £ £ 


x . r->exp<-iaj m t + iN d>) + cc 
mn - m n T 


( 2 ) 


where cc or * denotes the complex conjugate. We consider here the interac- 
tions among these frequency and azimuthal components along with their interac- 
tions with the mean flow and the background turbulence. These interactions 
are viewed within the framework of the energy exchanges among the different 
scales of motion. To derive the governing equations for each flow component, 
the usual phase-averaging technique (refs. 19 and 20) is extended to apply to 
a given azimuthal number in addition to the usual definition with respect to 
the frequency. Thus, 
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The mean flow momentum equation is derived by substituting the decomposi- 
tion of equation (1) into the full, unsteady Navier-Stokes equations for incom- 
pressible flow in cylindrical coordinates and by taking the time average. The 
mean flow energy equation is obtained by multiplying each component of the mean 
flow momentum equation by the corresponding velocity and adding. The result- 
ing mean flow kinetic energy equation is 
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where 
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The kinetic energy equation of the turbulent motion is derived as fol- 
lows. The time-averaged momentum equation is subtracted from the full momen- 
tum equation. Each component of the momentum equation is multiplied by the 
corresponding velocity component u ! , and the three components are added and 
time averaged. Triple correlations'that contain an odd number of random compo- 
nents are zero; those containing an even number of random components are not 
zero. The resulting kinetic energy equation of the turbulence is 


3 



IK - ll-V W 2 + W w'v' - Vw'w' + Wi 


j 3x.. r 


Ww' v 1 


9Ul / \ 3Ul g / ^ \ 

- u i u j 5^ * C- u i u j - u i u j) 3^ - 5*J ( u i <p ' * q) ♦ u j<0 


Re 


,2 m1 1 


u i 7 u i - 7 y 1 - 2V 


. awl _ JL w .2 _ 2w , 3Y1 

rt. oi” 


9<(> 


9<J> 


( 4 ) 


where 
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The energy equation for a wave component of frequency u m and azimuthal 
number N n Is obtained as follows. The full momentum equation Is first phase 
averaged with respect to the frequency under consideration, and the time- 
averaged equation Is subtracted from the phase-averaged one. The resulting 
equation contains all the azimuthal components of frequency w m . To separate 
the azimuthal component, we multiply the w-momentum equation by exp(-1N n <j>) and 
Integrate It over <j>. This produces a momentum equation for the mn-frequency- 
azlmuthal component. This equation Is multiplied by the corresponding velocity 
tfl.mn to yield the kinetic energy equation for the mn-wave component: 
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The dynamic equation for the wave-induced stresses rn mn Is derived as 
follows. The uj-momentum equation Is multiplied by u i , ana vice versa, to 
form the equations for the Instantaneous Reynolds stresses. By taking the 
difference between the phase average and time average of the instantaneous 
Reynolds stresses, the following equation for rjj Is obtained: 
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Integral Energy Equations 


The Integral energy equation of the mean flow is obtained by applying the 
boundary-layer-type approximations to the mean quantities. These approxima- 
tions imply that d( )/9x << 9( )/9r, V << U, and that the mean flow is 
axi symmetric with W = 0. These approximations are applied to equation (3), 
viscous terms are neglected, and the equation is integrated over r and 4>. 
The resulting mean flow energy equation is given by 
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As it will be shown later on, uv, in general, contains azimuthal terms 
that can distort the axi symmetry of the mean flow. These azimuthal terms van- 
ish upon integration over <j>. Therefore the integration over 4> was necessary 
to be consistent with the assumed axi symmetry of the mean flow. Equation (7a) 
therefore represents the mean flow energy in a transverse section of the jet. 

Equation (7a) can be written as 
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where 0 is the momentum thickness and MW and MT are the first and second 
integrals in the right side of equation (7a), respectively. Equation (7) 
states that the growth of the momentum thickness (the drain of the mean flow 
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energy) is governed by the mean flow production of the wave components MW and 
by the mean flow production of the turbulence MT. 

The integral energy equation for the random turbulence is obtained by 
applying boundary-layer- type approximations to the mean quantities and by han- 
dling the viscous dissipation terms in the usual manner. After integrating 
over r and <J>, the turbulence kinetic energy equation reduces to 
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This equation can be written as 
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Equations (8a) and (8b) state that the development of the turbulence energy is 
governed by its production by the mean flow MT, the energy transfer from the 
wave WT, and the turbulence dissipation dis. 

The integral energy equation for the mn-wave component is obtained from 
the full energy equation ( 5 ) in the same manner: 
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This system of equations can be written In the form 
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which states that the development of the wave energy is governed by the mean 
flow production of this particular wave MW mn , by the energy transfer between 
this wave and the turbulence WT mn , and by the interactions between this wave 
and other existing waves WW mn . 


Mode Decomposition and Interactions 

The integral energy equations (eqs. (7) to (9)) are the basis for the 
present analysis. In these equations, terms representing wave interactions 
are given in their most general forms, which are applicable to any number of 
frequency and azimuthal components. The decomposition of these terms for a 
given number of wave components yields 

u. = t' u, m „exp(-1u m t + iN„<j>) + cc 
i L-d i mn r m n T 
m,n 

In the mean flow energy equation (7), the production of all waves is given by 


The waves' Reynolds stress 

uv = > (u v* + 

Z-f V mn mn 

m,n m,k,!l 

The first summation is over the product uv* of the same wave. These terms 
are axisymmetric irrespective of their azimuthal number. The second summation 
over the product uv* is produced by the waves of the same frequency but dif- 
ferent azimuthal numbers. These terms are not axisymmetric and therefore 
destroy the symmetry of the mean flow. Note that these terms exist only if 
there is more than one azimuthal component of the same frequency. Thus, this 
is consistent with Cohen and Wygnanski's (ref. 21) conclusion that single mode 
excitation cannot destroy the axi symmetry of the mean flow. The azimuthal 
terms appearing in equation (11) vanish when integrated over <j>. Therefore 
these terms redistribute the mean flow energy in the azimuthal direction but 
do not contribute to the total energy integrated over 4>. Since the mean flow 
energy from equation (10) is integrated over <J>, these nonaxi symmetri c terms 
vanish and the mean flow production of the wave can be written as 
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Thus, the mean flow productions of the waves are given by the linear superposi- 
tion of the Individual mean flow productions of each wave. 

In the turbulence energy equation (8), the energy exchange between the 
waves and the turbulence WT Is given by 
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where the wave-induced stresses r . ^ ■ <u;u'.> - u ; u ' . Now let us decompose r 
In a manner similar to u: that Is, ^ ^ ^ 
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Thus, ru are generally given by 
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The first summation Is over the product produced by the same wave. The second 
summation Is over the product produced by two waves of the same frequency and 
different azimuthal numbers. This second summation vanishes upon Integration 
over <j>. Thus, with the energy defined as Integrated over <j>, we can write 
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Therefore, the turbulence energy exchange with the waves Is given by the sum of 
the energy exchanges between the turbulence and each Individual wave. 
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Next, we consider the wave-wave interactions appearing in equation (9): 
namely, 
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Because of the time averaging and the azimuthal integration, these terms are 
generally zero unless w and N simultaneously satisfy one of the following 
sets of conditions: 


or 
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These are the same relations given by Cohen and Wygnanski (ref. 21). Unfortu- 
nately, this is not quite a restrictive condition because if one starts with 
only two frequency components, they can amplify other frequency components 
that satisfy the conditions of equation (18). If the process is repeated, an 
enormous number of frequency components can be amplified. The present analy- 
sis is applicable to an unlimited number of components. However, for clarity 
we present only the wave-wave interactions for a limited number of wave compo- 
nents. We consider six wave components: three axi symmetric frequency compo- 

nents and three first-helical frequency components. The frequencies of the 
helical components match those of the axisymmetric components. These initial 
wave components can interact with each other to generate an enormous number of 
other frequency and azimuthal components. However, only the interactions among 
these initial components are considered. Thus, we force other frequency and 
azimuthal components to be Identically zero. 


The frequencies are chosen to satisfy the harmonic relations: that is 
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with m = 1,2,3. Thus, each two consecutive frequencies are related to each 
other by subharmonic-fundamental relations. The lower frequency denotes that 
of the subharmonic, and the higher denotes that of the fundamental. 

The conditions given by equation (18) on the frequency and azimuthal 
number limit the number of nonzero wave-wave interactions. The wave-wave 
interactions can be classified into three groups: (1) interactions among the 

axisymmetric waves, (2) interactions among the helical waves, and (3) mixed 
interactions between the axisymmetric and the helical waves. The general forms 
of these interactions will be discussed. The interactions among the axisymmet- 
ric waves can be manipulated to be in the form of subharmonic-induced stresses 
multiplied by fundamental-induced strains. The Interactions between the first- 
helical and the axisymmetric components are formed by two subharmonic frequency 
components of different azimuthal numbers interacting with the first-helical 
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fundamental component. The interactions among the helical modes produce non- 
axisymmetric terms that vanish upon integrating over <j>. Thus, with the energy 
defined as that integrated over <{>, the integrated interactions among the 
helical modes are Identically zero. This is an interesting result since it 
indicates that, unlike the axi symmetric modes, two first-helical frequency com- 
ponents cannot generate other frequency components. The energy exchanges among 
the several scales of motion are represented schematically in figure 1. The 
wave-wave Interactions for the six components considered are given by the fol- 
lowing equations: 

For the axi symmetric w-component, 

””l0 = -Wio, 10,20 + ”21 ,11,10 
For the axi symmetric 2aj-component, 

””20 - ”10,10,20 " ”20,20,30 + W 31 ,21,20 
For the axi symmetric 4w-component, 

” w 30 - ”20,20,30 

For the first-helical u-component, 

””ll * -”10,11,21 ~ ”21 ,11 ,10 
For the first-helical 2w-component, 

””21 = ”10,11,21 " ”20,21,31 “ ”31,21,20 
For the first-helical 4«-component, 

WW 3 1 = ”20,21 ,31 
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Closure Models 

The integral energy equations are used to obtain a set of ordinary differ- 
ential equations representing the energy exchanges among the different scales 
of motion. Each kind of motion is assumed to be descrlbable with a limited 
number of shape parameters. The shape assumptions for the mean flow, waves, 
and turbulence required to perform the radial integrals in the energy equations 
are discussed next. 

Mean flow . - We describe the mean flow in terms of the momentum thick- 
ness 9. The mean flow profile is the two-stage hyperbolic tangent profile 
proposed by Michalke (ref. 22). 

For 9 < 0.08, 

I I for 0 < r < 1 - |- 

\ |l + tanh[^ (1 - r)] | for 

For 0 > 0.08, 

U(r ,0) = !{l + tanhj^Q (y - r)] } for 

This profile fits the experimental data (e.g., refs. 23 and 24) for the mean 
flow development between the jet exit and the fully developed region. 


1 - |"< r < • > (22) 


0 < r < -» 
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The mean flow is thus characterized by the momentum thickness 9 rather 
than by the axial distance x. Therefore functional dependency on x in the 
integrals involved in the energy equations will be replaced by 9. 

Random turbulence . - The radial distribution of the turbulent stresses is 
assumed to be given by a Gaussian profile in the form 

2 

ujuj = a^ E(x) G(9)exp(-n ) 

where 

2 (r - l ) 2 

n = 5 


This assumption is consistent with Lau's (ref. 25) findings based on correlat- 
ing a large set of axial velocity fluctuations of unexcited jets over an 
extended range of jet operating conditions. The term E(x) is the turbulence 
energy, and a^j and C are constants given by 

ai i = 1 a 22 = S 33 * 0.5 a 1 2 * 0.33 a 3 3 = a 23 = 0 and C * 20 



(24) 

This model is similar to that used by Mankbadi and Liu (ref. 6 ) and by Mankbadi 
(ref. 7), but it has been modified to account for a given turbulence level in 
the potential core. Within the potential core, the turbulence level is deter- 
mined by its initial value at the exit of the jet. A normalization function 
G(9) renders the fine-grained turbulence energy E(x) over a section of the 
jet. Thus, 



for 9 < 0.08 and T 0 = 0. 


For 9 > 0.08 and T = 0.5 u'2(1 - 9/0.08)2, u <2 j S the initial turbulence 

o oc oc 

intensity at the center of the jet exit. 

This turbulence model is based on quasi-steady assumptions in which data 
from steady-flow experiments are used to model the turbulence at unsteady con- 
ditions. However, E(x) is determined later on from the nonlinear interac- 
tions. The experimental data thus provide only the radial shape and the ratios 
between the stresses. The shape and ratios may vary under excitation at high 
frequencies. However, the interest here is within the jet column mode, the 
so-called preferred frequencies, which in general are lower than those of the 
shear layer instability frequencies. Therefore, for this range of frequencies 
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it is logical to assume that the ratios of the turbulent stresses are given by 
their corresponding values for the steady-state case. 

The viscous dissipation in the turbulence energy equation is related to 
the small-scale turbulence as 


e 



(26) 


where &z - 1.5. 

In the present system of variables S is replaced by 0 where $ = a 3 © 
and a 3 = 4.4. 

Wave components . - The momentum equation for each wave component mn is 
obtained by subtracting the time-averaged momentum equation from the mn-phase 
averaged momentum equation: 
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where 6ij is Kroneker's delta. 

If this equation is linearized and the terms containing r^j are 
neglected, it can be simplified to the linear stability equation: 
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The shape assumptions for the wave components follow those of previous 
works (refs. 6, 7, 26, and 27) in assuming that the Fourier coefficient u mn 
in equation (2) can be separated into an amplitude that is a function of x 
and a shape that is a function of r at a given location along the jet. In 
these previous works, the complex radial shape was obtained from the linear 
stability equation. In Mankbadi (ref. 7) the amplitude was obtained from the 
nonlinear analysis while its phase angle was obtained from the linear analysis. 

Experimental (refs. 17 and 28) and theoretical (refs. 15, 16, and 29) 
investigations have indicated that nonlinear wave interactions are highly 
dependent on the phase angles between the waves. Furthermore, the experimental 
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data of Mankbadi, Raman, and Rice (ref. 28) clearly indicate that the develop- 
ment of these phase angles is a highly nonlinear process. Therefore, previous 
shape assumptions for the wave components are modified to 
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The radial shapes u(r,0) are given as the eigenfunctions of the locally par- 
allel linear stability equation as before. However, now both the amplitudes 
A mn and the phase angles Y mn are determined from the nonlinear analysis. 


As pointed out before, the linearized form of the momentum equation for 
the mn-wave is simply the linear stability equation. Therefore, to a first 
approximation <j/(x) is given by exp(i 0 J x a(£) d£> where a is the complex 
wave number corresponding to the frequency w. The radial shapes u(r,0) and 
the corresponding complex wave numbers are obtained from the local linear sta- 
bility theory. The fact that the radial shapes of the waves can be obtained 
from the locally parallel stability theory was confirmed experimentally by 
Strange and Crighton (ref. 30) and Wygnanski and Petersen (ref. 11). 


The solution of the inviscid stability equation follows that given by 
Michalke (ref. 22). In the damped region, a complex contour is used for the 
numerical integration as in Mankbadi and Liu (ref. 6). The eigenfunctions are 
normalized in such a way as to render A mn the energy of the mn-wave in a 
section of the jet. The arguments of the complex eigenfunctions are normalized 
by taking the phase angle of u to be zero at the jet centerline. The phase 
angle Is therefore that of the streamwise wave velocity at the jet centerline. 


Wave-induced stresses . - When the full form of the r^j equations is 
linearized, it reduces to the same form as that of a singlemn-wave component. 
Therefore, it is modeled following Mankbadi and Liu (ref. 6) for the interac- 
tion between the single-frequency component and the random turbulence. Based 
on the linearized form, it can be shown that r^j takes the form 


i J ♦mn 


A ( x ) 
mn 


E<x>f 11,n.n exp [ ,< 1' !x) - V * ,N n+] 


+ CC 


(30) 


The closure model for the wave-induced stresses follows that of the mean turbu- 
lence quantities. Details of the model are given in Mankbadi and Liu (ref. 6). 
The modeled linearized equation for r^j is 
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where *^-l .445311/ 3r >0 Is the time scale for return to Isotropy, r^ ■ 
r xx + frr + f"<b<b > and ' and D Indicate differentiation with respect to r. 


Nonlinear Interaction Equations 

With the shape assumptions Introduced In equation (3), the energy equa- 
tions (eqs. (7) to (9>) can be Integrated across the jet to yield the following 
system of simultaneous ordinary differential equations. 

The mean flow momentum thickness 0(x) Is given by 
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The turbulence kinetic energy 


Is given by 
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In this system of equations, I represents an integral that is, in 
general, a function of the momentum thickness 0, frequency w, and the azi- 
muthal number N considered. In the mean flow equation (32), 1^ is the mean 
flow advection integral given by 
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The first term on the right side of equation (32) is the mean flow production 
of the turbulence, and the second term is the production of each Individual 
wave. The mean flow production integrals of the random turbulence and the wave 
component, I^j and Imh> respectively, are given by 
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where Re Indicates the real part. 

In the turbulence energy equation (33), Ija is the integral of the turbu- 
lence advection by the mean flow given by 


I MT (e) 



U exp(-n )r dr 


(39) 


The first term on the right side of equation (33) is the turbulence produced 
by the mean flow. The second term is the turbulence energy exchanged with the 
waves, and the last term is the viscous dissipation. 
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The Integrals Ihj and 


are given by 
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In the energy equation for the mn-wave (eq. (34)), 1^ is the Integral 
for the wave advection by the mean flow: 
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The first term on the right side of equation (34) is the wave production 
by the mean flow. The second term is the interaction with the turbulence. The 
last term in equation (34) is the interaction of the mn-wave with other waves. 
The term WW mn is actually the sum of the interactions of the mn-wave with 
other existing waves. It is given here only for waves with frequencies of to, 
2co, and 4to and with azimuthal numbers of 0 and 1. Thus, for these six 
waves WW mn takes the form 
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The phase-development equation (35) is obtained in a manner similar to that of 
Mankbadi (ref. 16). The Strouhal number S is defined as u mn d/(2irl) e ) , where 
d is the nozzle diameter; 1 ^ is given by 
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and WW^ n is identical to WW mn (eqs. (42) and (43)) except that it is the 

imaginary part, and therefore the cos( ) term is replaced by sin( ): that 

is, 
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This system of equations (eqs. (32) to (35)) is subject to the initial condi- 
tions at x = 0 : 0(0), E -^ ( 0 ) , and E mn ( 0 ), and to the initial phase angle 
with respect to a given reference. 


TWO-FREQUENCY EXCITATIONS 

Before we present the results for multiple frequency excitations, we 
first compare the predictions of the present formulation with corresponding 
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experimental data for the case of a turbulent round jet excited at two frequen- 
cies. Linder two-frequency excitations, control can be exercised over the ini- 
tial energy levels as well as over the initial phase difference between the 
two excitation components. 

The solution of the system of nonlinear equations depends on the speci- 
fied Initial levels. Therefore, to compare the present predictions to the 
experimental results, we obtain the initial conditions from the corresponding 
experiment. 

The calculated fundamental and subharmonic components at Strouhal numbers 
of 0.3 and 0.6 are shown in figure 2 for several initial phase angles. When 
only two waves are present, there is one initial phase-difference angle 3 
between them. The initial momentum thickness is 0.026, and E 0 = 0.0001. The 
initial energies of the fundamental (S = 0.6) and the subharmonic (S = 0.3) 
are taken such that the initial instability axial velocity components at the 
jet centerline are 1.2 and 0.6 percent of the jet exit velocity. These values 
were based on Bradley's data (ref. 13). Figure 2 shows that the fundamental 
is not as sensitive to the phase difference as the subharmonic is. Bradley 
measured the spectral amplitudes that were proportional to the energies 
(ref. 13). His measured data behave similar to the predictions of figure 2. 

The fundamental is less dependent than the subharmonic on the initial phase 
angle. Bradley's measured maximum subharmonic amplification occurs at 
3 = 180°, and the minimum subharmonic amplification occurs at 3=0° 

(ref. 13), which is the same as the predicted results of figure 2. 

Figure 3 compares the calculated centerline phase-averaged velocities with 
the corresponding data of Arbey and Ffowcs-Wi 1 1 lams (ref. 18). The Strouhal 
numbers are 0.3 and 0.6, and E 0 = 0.00001. The initial centerline velocity of 
the S = 0.3 component is 1.5 percent of the jet exit velocity; the initial 
centerline velocity of the S = 0.6 component is 0.38 that of the S = 0.3, as 
in the experiment. At S = 0.3, the calculated peaks occur further downstream 
than the data peaks (fig. 3(a)). However, the calculated peaks are at the same 
velocity level as the data peaks. Both the theory and the experiment indicate 
that the S = 0.3 component peak increases when 3 changes from 0° to 180°. 

The calculated phase-averaged velocities at S = 0.6 (fig. 3(b)) also have 
the same features, the same level of amplification, and the same dependency on 
the phase angle as the measured data. The measured component increases again 
after it decays, which is probably due to its Interaction with other frequency 
components - a mechanism which is not allowed here since frequency components 
other than S = 0.3 and S = 0.6 are forced to be zero. 

The dependency of the subharmonic amplification on the initial phase angle 
is shown in figure 4 for Strouhal numbers of 0.2 and 0.4. The calculations 
are shown in figure 4(a), and the corresponding measurements of Raman and Rice 
(ref. 12) are shown in figure 4(b). In both cases the velocities are normal- 
ized by the maximum attainable peak of the subharmonic. The initial levels are 
Uf 0 = 7 percent, u so = 0.5 percent, and E 0 = 0.005. At this high level of 
excitation the acoustic decay is dominant in the first half diameter (ref. 12). 
Therefore, Tam and Morris's (ref. 31) correction for the acoustic interaction 
is used in the first half diameter. Figure 4 shows that both the measurements 
and the theory predict a strong influence of the phase angle. The subharmonic 
peak can be reduced to about one-third of its maximum attainable value depend- 
ing on the initial phase difference between the two waves. 
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The calculated and measured momentum thicknesses are shown In figure 5 
for the two-frequency excited and the unexcited cases. The initial conditions 
are the same as those of figure 4. Figure 5 shows that two-frequency excita- 
tion can considerably increase the jet mixing. However, the enhanced momentum 
thickness is only weakly dependent on the phase angle. This indicates that the 
direct role of the subharmonic in controlling the mixing of turbulent jets is 
less pronounced than its role in advancing the mixing of laminar jets (refs. 5 
and 16). However, the subharmonic can still have a strong indirect role in the 
mixing process since It can control the background turbulence that, in turn, 
controls the mixing process. 

If both the fundamental and subharmonic initial energy levels are large, 
their respective energy drain from the mean flow is also large. Therefore, 
the wave-wave Interactions become less significant. The dependency of the wave 
amplification on the initial phase angles should be less pronounced. This phe- 
nomenon is demonstrated in figure 6 where both waves have the same high initial 
energy levels, 0f o * u so * 3 percent. The Strouhal numbers are 0.3 and 0.6. 

The present predictions of the subharmonic peak as a function of the initial 
phase difference are shown in figure 6(a), and the corresponding data of Raman 
and Rice (ref. 12) are shown in figure 6(b). The theory and observations both 
Indicate the same weak dependency on the phase angle. 

According to the present analysis, the wave-wave interactions are depend- 
ent on the energy levels. Therefore, increasing the initial energy level of 
the fundamental, while keeping that of the subharmonic small, should Increase 
the peak of the subharmonic. The calculated subharmonic peak as a function of 
the initial level of the fundamental Is shown in figure 7(a). The Strouhal 
numbers are 0.2 and 0.4, the initial velocity of the subharmonic is 1/15™ 
that of the fundamental, and the Initial phase angle is kept at 270°. The cor- 
responding data of Raman and Rice (ref. 12) are shown in figure 7(b). Both 
theory and experiment indicate that the peak of the subharmonic Increases as 
the fundamental's energy level Increases. The subharmonic reaches a satura- 
tion value of about 20 percent. 


MULTI FREQUENCY EXCITATION 

The data of Moore (ref. 32) and Raman, Rice, and Mankbadl (ref. 33) indi- 
cate that within the Strouhal number range 0.2 to 1.0, a forced-wave component 
amplifies along the jet and can influence the growth rate of the jet. There- 
fore, in the present study three frequencies, corresponding to Strouhal num- 
bers of 0.2, 0.4, and 0.8, are considered. These Strouhal numbers were chosen 
to cover jet excitability in the experimentally indicated range referred to as 
the "jet column mode" (ref. 9). The next two sections present the results for 
multifrequency excitation in the axisymmetric mode. In the first section 
results for interactions with equal initial levels of all waves are presented. 
In the succeeding section, results for variable ratios of initial energy levels 
are presented. 


Multifrequency Forcing at Equal Energy Levels 

In this section, the Initial energy levels of all waves are equal, and the 
initial momentum thickness is taken to be 0.012. Since the focus here is on 
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highly turbulent jets, the Initial turbulence energy level is taken to be rela- 
tively high Et<0) = 0.005. With three interacting waves there are two initial 
phase differences. The phase difference between S = 0.2 and S = 0.4 is 
3i 2 and that between S = 0.4 and S * 0.8 is 323 • 

Figure 8 shows the energy peak of each wave component normalized by its 
initial value as a function of the phase differences. The initial energy 
level (which is equal for the three waves) is varied from 10 -8 to 5xl0 -8 in 
figure 8(a) to (d). This range of initial energy levels corresponds to ini- 
tial wave velocities of roughly 0.001 to 0.07 U e . A comparison of part (a) of 
figure 8 with parts (b) to (d) indicates that the dependency on the phase 
angle Increases as the Initial energy level increases. This trend is expected 
since the wave-wave interaction is a nonlinear process. At low levels 
(fig. 8(a)), each wave grows almost independently of the other waves. At 
higher energy levels (fig. 8(c)), the energy peak can be doubled depending on 
the phase angle, indicating strong nonlinear interactions. 

A comparison of parts (a) to (d) of figure 8 also indicates that the maxi- 
mum attainable normalized peak decreases as the initial energy levels increase. 
This fact is demonstrated in figure 9, which shows the maximum attainable peak 
(depending on the combination of phase differences) as a function of the ini- 
tial energy level. Figure 9 indicates that at low initial energy levels the 
peak varies in accordance with the linear theory. However, the rate of growth 
decreases with higher energy levels, and saturation conditions are reached 
around energy levels of 0.005 to 0.01. Further increases in the initial energy 
level damp the forced component. 

The maximum attainable (depending on the phase angles) momentum thickness 
along the jet is shown in figure 10 for several initial energy levels. This 
figure shows that the higher the initial energy level, the greater the mixing 
rate enhancement due to excitation. Since we are Interested in enhancing the 
mixing and since the saturation is limited as shown in figure 9, E 0 = 0.005 
seems to be the optimum value. Therefore, multifrequency forcing at 
E 0 = 0.005 is examined further in figures 11 to 21. 

To examine the role of multifrequency forcing on jet mixing, we use three 
phase-difference combinations but keep the other initial conditions fixed 
^igs. 11 to 16). In figures 11 and 12 the phase differences U 12 and 323 
are 0° and 90°, respectively. Figure 11 shows the energies of the three fre- 
quency components and indicates that S = 0.4 is the most amplified component 
along the jet. The mean flow productions of the wave components and of the 
turbulence are shown in figure 12. The sum of these productions determines 
the growth rate, and hence the mixing rate, of the jet. The mean flow produc- 
tion of the turbulence for the unexcited case is also shown in figure 12. Fig- 
ure 12 shows that the production of the three waves is important in the initial 
region of the jet. But further downstream beyond x = d, the production of the 
turbulence dominates. 

In figures 13 and 14 the phase angles 3i 2 and 323 are 90° and 180°, 
respectively. The energies of the waves are shown in figure 13, which indi- 
cates that the S = 0.4 component is less amplified than it is in figure 11. 

The mean flow productions are shown in figure 14, which also indicates that the 
turbulence production is the dominant mechanism beyond x = d. The case of 
3l 2 ■ 270° and 323 = 180 ° is shown in figures 15 and 16. Figure 15 shows 
that for this combination of phase differences S = 0.2 is the most amplified 
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component along the jet. The corresponding mean flow productions are shown in 
figure 16. This figure also shows that turbulence production is the dominant 
mechanism. A comparison of the turbulence produced by the mean flow for the 
excited and unexcited cases (fig. 16) shows that multifrequency excitation can 
considerably increase the random turbulence produced by the mean flow. 

The momentum thickness along the jet is shown in figure 17 for the three 
pairs of phase-difference combinations that were presented in figures 11 to 16. 
Figure 17 shows considerable enhancement in the momentum thickness for the 
three cases considered. In the case of £12 = 90° and 023 * 0°, the momen- 
tum thickness around x = 4d is greater than for 0]2 * 90° and 023 = 180 ° 
because the S = 0.4 component is amplified more in the former case than in the 
latter case (fig. 11 versus fig. 13). Figure 17 also shows that the momentum 
thickness is most enhanced when 2 * 270° and 023 =180°. For this case 
S * 0.2 is the most amplified component (fig. 15). This low Strouhal number 
wave has the longest streamwise lifespan, and therefore it is the most effec- 
tive in enhancing mixing at downstream locations. 

Figure 18 shows the turbulence energy along the jet for the phase- 
difference combinations in figure 17. The turbulence energy for the unexcited 
case is also shown on figure 18. The turbulence is enhanced because of the 
forcing. Note also that forcing at different phase-difference combinations 
has the same effect on the turbulence as it does on the momentum thickness 
(fig. 17). 

In laminar jets, where turbulence is negligible, the growth rate is dic- 
tated by only the production of the waves. Therefore subharmonic amplification 
or damping would have a pronounced effect on the jet growth rate (refs. 5 and 
16). For turbulence jets where the growth rate is governed by the production 
of the turbulence and the waves, the situation is different. As indicated by 
figures 12, 14, and 16, the turbulence production is the dominant mechanism in 
controlling the growth rate. The direct role of wave amplification in the 
growth rate is not pronounced for the turbulent jets as it is for the laminar 
jets. The forced waves, however, enhance the random turbulence and therefore 
still play an indirect role in the mixing process. 

The turbulence enhancement of the random turbulence that was shown in 
figure 18 is evident in the early experiments of Binder and Favre-Marinet 
(ref. 3). The mechanism of turbulence enhancement under multifrequency excita- 
tion is examined here. The energy pumped into the turbulence by the waves and 
the mean flow are shown In figures 19 to 21 for three different pairs of phase- 
difference combinations. These figures show that turbulence enhancement under 
excitation is a result of a direct and an indirect mechanism. In the direct 
mechanism, the imposed waves decay by dissipating their energy into the back- 
ground turbulence (the second term in the turbulence energy equation, 
eq. (33)). In the second indirect mechanism, the wave acts as a catalyst that 
enhances the mean flow production of the turbulence. Note that in figures 19 
to 21, the mean flow production of the turbulence for the excited case is much 
higher than that in the unexcited case. The mean flow production of the turbu- 
lence is proportional to the local level of the turbulence energy (the first 
term in the turbulence energy equation (eq. (33))). Therefore, wave-produced 
enhancement in E* further increases the mean flow production of the turbu- 
lence EtlMT- Thus, the mean flow production of the turbulence is increased 
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by excitation. In figure 21, the S = 0.2 component, which has a long stream- 
wise lifespan, is amplified. This wave keeps pumping energy to the turbulence 
for a long distance along the jet and causes the greatest enhancement in the 
mean flow production of the turbulence. 


Forcing at Variable Initial Energy Ratios 

In the previous section we considered the case of forcing at equal energy 
levels. In the next section, we consider the case when one of the waves has an 

initial energy level that is one order of magnitude lower than the other two. 

In the succeeding section, we consider the case when one of the waves has an 
initial energy level that is one order of magnitude higher than the other two. 

One wave initially weaker than the other waves . - The peak of the three 
waves as a function of the initial phase angles are shown in figure 22. In 
figure 22(a) the S = 0.8 wave is initially an order of magnitude smaller than 

the other two. In figure 22(b), the S * 0.4 wave is the weak one, while in 

figure 22(c), the S = 0.2 wave is the weak one. In figure 22(a), the S = 0.8 
component is negligible with respect to the other two waves. Therefore the 
S = 0.2 and S = 0.4 waves behave as in the case of two-frequency excitation. 
When S = 0.4 is the weak wave, figure 22(b), it still amplifies considera- 
bly. The S = 0.4 wave has a Strouhal number near that of the most amplified 
wave. Furthermore, this wave interacts directly with the other two waves. 
Considerable wave Interactions occur that can amplify the S » 0.4 wave to the 
same level as that of the other two waves. When S = 0.2 is the weak wave, 
the S = 0.2 wave is not amplified and the other two waves, S = 0.4 and 
S = 0.8, behave as if the third wave were not present at all (fig. 22(c)). 

The turbulence energy along the jet is shown in figure 23 for the initial 
energy ratios of figure 22. Note that the figure Indicates that the S = 0.2 
component must be present to have a considerable effect on the turbulence 
because of the long streamwise lifespan of the S = 0.2 component, as indicated 
before. The figure also shows a greater effect on the turbulence when the two 
high-level waves are S = 0.8 and 0.4, instead of S = 0.2 and 0.4. As indi- 
cated before, when the S = 0.8 and S = 0.2 waves are high, they excite the 
S = 0.4 component and thus act effectively as three-frequency excitation. 

This mechanism is not present when the two high-amplitude waves are S = 0.2 
and 0.4. 

The momentum thickness along the jet is shown in figure 24 for the same 
initial values as in figure 23. Figure 24 shows that the effect of excitation 
on the momentum thickness follows the same pattern as the effect of excitation 
on the turbulence. This indicates again that it is the turbulence that con- 
trols the mixing enhancement. 

One wave initially higher than the other waves . - One wave was given a 
high initial energy level of E 0 = 0.005 while the other waves were kept an 
order of magnitude lower (i.e., E 0 = 0.0005). Although this is a case of mul- 
tifrequency interactions, it can represent the experimental situation when only 
one frequency is forced. In an unexcited jet, there are several stability com- 
ponents that are naturally present. In an excited jet, the excitation wave can 
interact with other naturally present waves at other frequencies and amplify 
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them (e.g., ref. 9). The low initial energy level of E 0 = 0.0005 can repre- 
sent that of the naturally present waves. 

Figure 25 shows the energy peak for each frequency component as a function 
of the phase-difference angles. In figure 25(a) the S = 0.8 component is the 
high-amplitude wave. It produces some amplification of the S = 0.4 component, 
which in turn can enhance the S = 0.2 component. In figure 25(b), the S = 0.4 
component is the high-amplitude wave. It can enhance the S = 0.2 component if 
the phase difference between them is 270°, but it has a negligible effect on 
the S = 0.8 component. In figure 25(c), the S = 0.2 component is the high- 
amplitude wave. The S = 0.4 component interacts with both the S = 0.2 compo- 
nent and the slightly amplified S = 0.8 component. 

The turbulence energy is shown in figure 26. The initial phase angles 
that produced the greatest effect on the turbulence were chosen for a given 
high-amplitude excitation at a given frequency. At x < d, the S * 0.8 compo- 
nent produces the highest effect on the turbulence. The S = 0.8 component 
amplifies very close to the jet exit and therefore is effective in pumping 
energy from the mean flow to the turbulence in this region. If one considers 
the region x > d, the S * 0.2 component is the most effective in pumping 
energy from the mean flow to the turbulence. The S = 0.2 component peaks fur- 
ther downstream and has the longest streamwise lifespan. Therefore, the high- 
frequency components are effective in enhancing the turbulence in the initial 
region of the jet while the low-frequency components are effective further 
downstream. The momentum thickness for the cases in figure 26 are shown in 
figure 27. Figure 27 shows that the effect of frequency on the momentum thick- 
ness follows the effect of the frequency on the turbulence. This again indi- 
cates that the turbulence is the effective mechanism for turbulence enhancement 
of excited turbulent jets. 

Figures 28 to 33 show the energy of the frequency components and the cor- 
responding turbulence productions for the three cases presented in figures 26 
and 27. In figure 28, S = 0.8 is Initially the highest component; it ampli- 
fies in the initial region of the jet. Figure 29 shows that the S = 0.8 wave 
component pumps considerable energy into the turbulence in the initial region 
of the jet. As a result of this energy transfer, the mean flow production of 
the turbulence is also enhanced. When the high energy is in the S = 0.4 compo- 
nent (figs. 30 and 31), the behavior is similar to that when S = 0.8 is the 
high component. However, the high S = 0.4 component peaks further downstream. 
As a result, the enhancement of the mean flow production of the turbulence in 
figure 31 is delayed in comparison with the corresponding enhancement in fig- 
ure 29. Note also that when S = 0.4 is the high-amplitude component, it 
dampens the S = 0.8 component (fig. 30). Therefore the S = 0.8 component plays 
a negligible role in pumping energy from the mean flow to the turbulence. When 
S = 0.2 is the high-amplitude component, the other frequencies are amplified 
slightly (fig. 32). In figure 33 all waves pump energy to the turbulence. 
Furthermore, because of its long streamwise lifespan, the S = 0.2 component 
keeps pumping energy to the turbulence along the jet. As a result of this con- 
tinuous increase in the turbulence energy, the mean flow production of the tur- 
bulence is also considerably enhanced. Apparently, if the S = 0.2 component 
can be amplified through its Interaction with other frequency components, it 
can have a considerable effect on the turbulence. 
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CONCLUSIONS 


The interactions among several stability waves in a turbulent round jet 
were examined based on the integrated energy in a section of the jet for the 
different scales of motion. The present formulations indicate that if only 
axisymmetric wave components are present, two frequency components with high 
amplitudes can interact with each other as well as with other background- 
frequency components and can generate an enormous number of other frequency 
components in the axisymmetric mode. If the frequency components present are 
only in the first-helical mode, two frequency components in the first-helical 
mode cannot interact directly with each other or with other first-helical fre- 
quency components. However, if both the axisymmetric and the first-helical 
component are simultaneously present with frequency components related by 
fundamental-subharmonic relations, mixed interactions between the axisymmetric 
and the helical mode occur. These mixed interactions amplify other frequencies 
in both the axisymmetric and the helical modes. 

The present predictions for two-frequency axisymmetric excitation produce 
results consistent with those from several published experiments. When a jet 
is excited at fundamental and subharmonic frequencies with the fundamental 
level higher than that of the subharmonic, the amplification of the subharmonic 
greatly depends on the initial phase difference between the two waves. This 
dependency on the initial phase angle becomes less pronounced if both compo- 
nents have high initial energy levels. Two-frequency excitation enhances the 
momentum thickness in accordance with observations. 

If a jet is forced with multifrequency components of equal initial energy 
levels, the nonlinear process depends on the initial level. At low energy lev- 
els, nonlinear interactions are negligible, and each wave behaves as if the 
others were not present. Increasing the initial energy level equally for all 
waves results in considerable nonlinear Interactions. At much higher initial 
levels, a saturation condition is reached in which the waves can no longer 
amplify each other. Increasing the initial energy level increases the momentum 
thickness of the jet until the saturation condition is reached. 

Examining the energy exchanges at high forcing levels of equal magnitudes 
for all the waves indicates several Interesting features. The enhancement in 
the momentum thickness is dominated by the mean flow turbulence production, 
rather than by the fundamental-subharmonic interactions as in the laminar case. 
Multifrequency forcing of the jet enhances the turbulence through two mechan- 
isms. In the direct mechanism, the imposed waves decay by dissipating their 
energy into the background turbulence. In the second indirect mechanism, the 
wave acts as a catalyst that enhances the mean flow production of the turbu- 
lence. The turbulence produced by the mean flow is proportional to the local 
level of the turbulence energy. Therefore, wave-produced enhancement in the 
turbulence energy further increases the turbulence produced by the mean flow. 
Thus, excitation enhances the mean flow production of the turbulence. 

The role of the frequency on the nonlinear process was examined by 
varying the ratios of initial energy levels between the frequency components. 
Two-frequency excitation at high and low frequencies can result in considerable 
amplification of the intermediate frequencies. Thus, proper two-frequency 
excitation can be as effective as three-frequency excitation. 
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By taking the initial energy level of one of the frequency components to 
be much larger than the others, we demonstrated that high frequencies are 
effective in pumping energy from the mean flow to the turbulence. Low fre- 
quencies have longer streamwise lifespans and therefore are more effective in 
pumping energy from the mean flow to the turbulence at locations further 
downstream. 


APPENDIX - SYMBOLS 


A 

amplitude 

a i j > C 

constants 

cc 

complex conjugate 

D 

differentiated with respect to r 

d 

nozzle diameter 

dis 

turbulence di ssipation 

E(x) 

energy 

E 0 

initial energy level 

E P 

energy peak 

Et 

turbulence energy 

G<e) 

normalization function 

i 

integral that is a function of 0, co, and N 

IMA 

mean flow advection integral 

Imt 

mean flow production integral of the random turbulence 

Imw 

mean flow production integral of the wave component 

J TA 

turbulence advection integral 

J WA 

wave advection integral 

9m 

imaginary 

9 . i ,j,k,St,m,n 

variables 

MT 

mean flow production of turbulence 

MW 

mean flow production of wave components 

N 

azimuthal number 
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p 

0 

q 

R 

Re 

Re 

r 

Hj.mn 

S 

T-l 

T 

o 

t 

U 

Oi Cx) 

u p 

u 

u! (x,r) 
u 

Uf 

M i ( X , t ) 

V 

V 

W 

WT 

WW 

w 

X 


pressure 

kinetic energy of the wave 
kinetic energy of the turbulence 
nozzle radius 
Reynolds number 
real part 

radical coordinate 
wave-induced stresses 
Strouhal number 

time scale for return to Isotropy 

0.5 u'2(l - 6/0.08) 2 
oc 

time 

velocity 

time-averaged motion 
peak velocity 
streamwise velocity 

background, fine-scale random turbulence 
periodic component 

periodic component at fundamental frequency 
periodic, organized, large-scale wave-shaped motion 
mean radial velocity 
radial velocity 
azimuthal velocity 

energy transfer between the waves and the turbulence 
interaction between one wave and other existing waves 
azimuthal velocity (angular velocity) 
axial distance 
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a 

13 

&12 

323 

6 

6 U 

e 

0 

P 

a 

<t> 

Y 

a) 

Subscripts: 

e 

f 

i ,j,k,a,m,n 

MA 

MT 

o 

oc 

s 

sp 

TA 

t 

WA 

x 


complex wave number corresponding to frequency u> 

initial phase-difference angle 

initial phase difference between S] and S2 

initial phase difference between S2 and S3 

displacement thickness 

Kroneker's delta 

viscous dissipation 

momentum thickness 

fluid density 

argument of I 

azimuthal angle 

phase angle 

frequency 

jet exit 

fundamental frequency 
variables 

mean flow advection 
mean turbulence 
initial conditions 
exit, centerline 
subharmoni c 
subharmonic peak 
turbulence advection 
turbulence (p. 28 ) 
wave advection 
axial coordinate 


31 



Superscripts: 


usual time average of the flow quantity 
wave 1 Ike 


complex conjugate 
eigenfunction 

differentiated with respect to r 
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(a) STROUHAL NUMBER, S, 0.3. <b) STROUHAL NUMBER, S, 0.6. 

FIGURE 3. - COMPARISON BETWEEN CALCULATED PHASE-AVERAGED CENTERLINE VELOCITIES 
AND ARBEY AND FFOWCS-WILLIAMS' DATA (REF. 16). 




(b) EXPERIMENTAL DATA OF RAMAN AND 
RICE (REF. 10). 


FIGURE 4. - DEPENDENCE OF SUBHARMONIC AMPLIFICATION ON THE INITIAL PHASE 
ANGLE, 0, AT STROUHAL NUMBER PAIR 0.4 AND 0.2. 
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(3) CALCULATIONS. (b) EXPERIMENTAL DATA OF MANKBADI, 

RAMAN, AND RICE (REF. 26). 


FIGURE 5. - DEVELOPMENT OF MOMENTUM THICKNESS UNDER TWO-FREQUENCY EXCITA- 
TION AND FOR THE UNEXCITED CASE. 




INITIAL PHASE ANGLE, p, DEG 

(3) THEORY. (b) EXPERIMENTAL DATA OF RAMAN AND 

RICE (REF. 10). 

FIGURE 6. - DEPENDENCY OF SUBHARMONIC PEAK ON THE PHASE ANGLE AT HIGH LEVELS OF 
BOTH FUNDAMENTAL AND SUBHARHONIC. 




FORCING VELOCITY, U f , PERCENT 

(a) THEORY. (b) EXPERIMENTAL DATA OF RAMAN AND 

RICE (REF. 10). 

FIGURE 7. - DEPENDENCY OF SUBHARMONIC PEAK ON THE FORCING LEVEL AT STROUHAL 
NUMBERS OF 0.2 AND 0.4. 


37 



NORMALIZED ENERGY PEAK, E n /E 



(b) INITIAL ENERGY LEVEL, E 0 = 10“ 4 . 



(d) INITIAL ENERGY LEVEL, E 0 * 0.5x10' 2 . 

FIGURE 8. - EFFECT OF PHASE DIFFERENCES, fl 12 AND p 2y ON PEAK 
OF EACH FREQUENCY COMPONENT (THREE STROUHAL NUMBERS, S). 
EQUAL INITIAL ENERGY LEVELS. 



LOG OF INITIAL ENERGY LEVEL, LOG (E p ) 

FIGURE 9. - MAXIMUM ATTAINABLE PEAK AS FUNCTION OF 
THE INITIAL ENERGY LEVEL. 



0 2 9 6 8 10 


x/d 

FIGURE 10. - EFFECT OF INITIAL ENERGY LEVEL ON THE DE- 
VELOPMENT OF THE MOMENTUM THICKNESS. 
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NORMALIZED ENERGY, NORMALIZED ENERGY, E/E 

E/E n 



FIGURE 11. - DEVELOPMENT OF THE ENERGY OF EACH FREQUENCY FIGURE 12. - MEAN FLOW ENERGY DRAIN TO THE WAVES (THREE 

COMPONENT (THREE STROUHAL NUMBERS, S) ALONG THE JET AT STROUHAL NUMBERS, S) AND TURBULENCE AT EQUAL INITIAL 

EQUAL INITIAL ENERGY LEVELS, E Q = 0.005. INITIAL PHASE ENERGY LEVELS, E Q = 0.005. INITIAL PHASE DIFFERENCES: 

DIFFERENCES: 0 12 = 90 u ; 0 23 = 0°. P, 2 = 90°; P 2J = 0°. 



FIGURE 13. - DEVELOPMENT OF THE ENERGY OF EACH FREQUENCY 
COMPONENT (THREE STROUHAL NUMBERS, S) ALONG THE JET AT 
EQUAL INITIAL ENERGY LEVELS, E Q = 0.005. INITIAL PHASE 
DIFFERENCES: P 12 = 90°; ^ = 180°. 


FIGURE 14. - MEAN FLOW ENERGY DRAIN TO THE WAVES (THREE 
STROUHAL NUMBERS, S) AND TURBULENCE AT EQUAL INITIAL 
ENERGY LEVELS, E Q = 0.005. INITIAL PHASE DIFFERENCES: 

P 12 = 90°; P 23 = 180°. 



0 2 4 6 


FIGURE 15. - DEVELOPMENT OF THE ENERGY OF EACH FREQUENCY 
COMPONENT (THREE STROUHAL NUMBERS, S) ALONG THE JET AT 
EQUAL INITIAL LEVELS, E Q = 0.005. INITIAL PHASE DIFFER- 
ENCES: P 12 = 270°; P 23 = 180°. 
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TURBULENCE LEVEL. E, PRODUCTION 



FIGURE 16. - MEAN FLOW ENERGY DRAIN TO THE WAVES (THREE FIGURE 17. - EFFECT OF INITIAL PHASE DIFFERENCES 

STROUHAL NUMBERS, S) AND TURBULENCE AT EQUAL INITIAL ( P 12 AND 3 23 > ON THE MOMENTUM THICKNESS. EQUAL 

ENERGY LEVELS, E 0 = 0.005. INITIAL PHASE DIFFERENCE, INITIAL ENERGY LEVELS, E = 0.005. 

P 12 = 270°. 


10x10‘ 3 


FIGURE 18. - EFFECT OF INITIAL PHASE DIFFERENCES, & 12 AND 
0 23 , ON THE TURBULENCE ENERGY. EQUAL INITIAL ENERGY 
LEVELS, E 0 = 0.005. 


FIGURE 19. - TURBULENCE PRODUCTION BY THE MEAN FLOW AND 
WAVES. EQUAL INITIAL ENERGY LEVELS, E Q = 0.005. INITIAL 
PHASE DIFFERENCES: 0 12 = 90°; & 23 = 0°. 
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MEAN FLOW 


o 



10x1 0" 3 



x/d 


FIGURE 20. - TURBULENCE PRODUCTION BY THE MEAN FLOW AND 
WAVES (THREE STROUHAL NUMBERS, S). EQUAL INITIAL ENERGY 
LEVELS, E 0 = 0.005. INITIAL PHASE DIFFERENCES: P 12 = 90°: 

P 23 = 180°. 


14X10* 3 



x/d 


FIGURE 21. - TURBULENCE PRODUCTION BY THE MEAN FLOW AND 
WAVES (THREE STROUHAL NUMBERS, S). EQUAL INITIAL ENERGY 
LEVELS, E 0 = 0.005. INITIAL PHASE DIFFERENCES: fy 2 = 
270° ; P 23 =180°. 
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ENERGY PEAK 




(3 ) INITIAL ENERGY LEVELS: E 10 = 0.005; E 20 




FIGURE 23. - TURBULENCE ENHANCEMENT UNDER VARIABLE RATIOS 
OF INITIAL ENERGY LEVELS. CASE A: INITIAL PHASE DIFFER- 

ENCES, p 12 = 270°; P 2 3 =0 ° ; ^ITIAL ENERGY LEVELS, E 10 * 
0.005; E 20 = 0.005; E 30 = 0.0005. CASE B: 0 12 = 270°; 

P 23 =0°; Eio = 0 - 005; E 2( r 0.0005; E 3Q = 0.005. CASE C: 

P 12 = 270°; p 23 = 270°; E 1q = 0.0005; E 2Q = 0.005; E 3() = 0.005 





(b) INITIAL ENERGY LEVELS: E 10 = 0.005; E 2Q = 0.0005; E 3Q = 0.005. 





(C) INITIAL ENERGY LEVELS: E 10 = 0.0005; E 20 = 0.005; E 30 = 0.005. 

FIGURE 22. - EFFECT OF INITIAL PHASE DIFFERENCES, 0 12 AND 0 23 , 

ON THE PEAK OF EACH FREQUENCY COMPONENT (THREE STROUHAL NUMBERS, 
S). NONEQUAL INITIAL ENERGY LEVELS. 



x/d 


FIGURE 24. - MOMENTUM THICKNESS ENHANCEMENT UNDER VARIABLE 
RATIOS OF INITIAL ENERGY LEVELS. CASE A: INITIAL PHASE 

DIFFERENCES. P 12 = 270°; P 23 = 0°; INITIAL ENERGY LEVELS, 
E 1q = 0.005; E 20 = 0.005; E 3Q = 0.0005. CASE B: P 12 = 

270°; p 23 = 0°; E 10 = 0.005; E 2Q = 0.0005; E 30 = 0.005. 
CASE C: P 12 = 270°; P 23 = 270°; E 10 - 0.0005; E 2Q = 0.005 

Ejq = 0.005. 
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ENERGY PEAK 




0 1234567 89 


FIGURE 26. - TURBULENCE ENHANCEMENT UNDER VARIABLE RATIOS 
OF INITIAL ENERGY LEVELS. CASE A: INITIAL PHASE DIF- 

FERENCES, P 12 = 270°; P 23 = 90 o ; INITIAL ENERGY LEVELS, 

E 10 = 0.0005; E 20 - 0.0005; E 3Q = 0.005. CASE B: 0 12 = 270°; 
P 23 = 270°; E 1q = 0.0005; E 2Q = 0.005; E 3() = 0.0005. CASE C: 
P 12 = 270°; P 2 3 = 180 ° ; E 1q = 0.005; E 20 = 0.0005; E 30 = 0.0005. 



0 123456789 


(C) INITIAL ENERGY LEVELS: E 10 = 0.005; E 20 = 0.0005; E 3Q = 0.0005. x/(J 

FIGURE 25. - EFFECT OF INITIAL PHASE DIFFERENCES, AND p 23 , ON FIGURE 27. - MOMENTUM THICKNESS ENHANCEMENT UNDER VARIABLE 

THE PEAK OF EACH FREQUENCY COMPONENT (THREE STROUHAL NUMBERS, S). RATIOS OF INITIAL ENERGY LEVELS. CASE A: INITIAL PHASE 

NONEQUAL INITIAL ENERGY LEVELS. DIFFERENCES, p 12 = 270°; P 23 = 90°; INITIAL ENERGY LEVELS, 

E 10 = 0.0005; E 20 = 0.0005; E 3Q = 0.005. CASE B: p 12 = 270°; 

fl 23 = 180°; E 10 = 0.0005; E 2Q = 0.005; E 3Q = 0.0005. CASE C: 
P 12 = 270°; P 23 = 270 o ; E 10 = 0.005; E 20 = 0.0005; E 3Q = 0.0005 
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FIGURE 28. - DEVELOPMENT OF THE FREQUENCY COMPONENTS AT FIGURE 29. - TURBULENCE PRODUCTIONS BY THE MEAN FLOW AND 

NONEQUAL INITIAL ENERGY LEVELS: E 10 = 0.0005; E 20 = WAVES AT NONEQUAL INITIAL ENERGY LEVELS: E 10 = 0.0005; 

0.0005; E 30 = 0.005. INITIAL PHASE DIFFERENCES: P 12 = E 2Q = 0.0005; E 50 = 0.005. INITIAL PHASE DIFFERENCES: 

270°; ^23 = 90°. 0 12 = 270°; 0 23 = 90 °- 


.015 

.010 


.005 


0 123456789 

x/d 

FIGURE 30. - DEVELOPMENT OF THE FREQUENCY COMPONENTS AT 
NONEQUAL INITIAL ENERGY LEVELS: E 10 = 0.0005; E 2Q = 

0.005; E™ = 0.0005. INITIAL PHASE DIFFERENCES: p 12 = 

270°; P 23 = 180°. 


0 123456789 

x/d 

FIGURE 31. - TURBULENCE PRODUCTIONS BY THE MEAN FLOW AND 
WAVES AT NONEQUAL INITIAL ENERGY LEVELS: E 10 = 0.0005; 

E 2Q = 0.005; E 30 = 0.0005. INITIAL PHASE DIFFERENCES: 
P 12 = 270°; P 23 = 180°. 





0 123456789 


x/d 

FIGURE 32. - DEVELOPMENT OF THE FREQUENCY COMPONENTS AT 
NONEQUAL INITIAL ENERGY LEVELS: E 10 * 0.005; E 2Q = 

0.0005; E™ = 0.0005. INITIAL PHASE DIFFERENCES: P 12 = 
270°; p 23 = 270°. 
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FIGURE 33. - TURBULENCE PRODUCTIONS BY THE MEAN FLOW AND 
WAVES AT NONEQUAL INITIAL ENERGY LEVELS: E 10 = 0.005; 

E 20 = 0.0005; E 30 = 0.0005. INITIAL PHASE DIFFERENCES: 
P , 2 - 270°: P 23 = 270°. 
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